Astronomy &: Astrophysics manuscript no. peddosc 


© ESO 2012 


August 15, 2012 





Oscillations of the Eddington Capture Sphere 

Maciej Wielgus 1 ' 2 , Adam Stahl 3 , Marek Abramowicz 3 ' 1,4 , and Wlodek Kluzniak 1 



(N 

o 

< 



1 Copernicus Astronomical Center, ul. Bartycka 18, PL-00-716 Warszawa, Poland 
e-mail: wlodek@camk.edu.pl 

2 Institute of Micromechanics and Photonics, ul. sw A. Boboli 8, PL-02-525, Warszawa, Poland 
e-mail: maciek.wielgus@gmail.com 

3 Physics Department, Gothenburg University, SE-412-96 Goteborg, Sweden 
e-mail: gusstaadOstudent .gu.se 

e-mail: marek . abramowicz@physics .gu.se 

4 Institute of Physics, Faculty of Philosophy and Science, Silesian University in Opava, Bezrucovo nam. 13, CZ- 746-01 
Opava, Czech Republic 

Received Aug 14, 2012; accepted ???? 



6 



ABSTRACT 

We present a toy model of mildly super-Eddington, optically thin accretion onto a compact star in the Schwarzschild 
metric, which predicts periodic variations of luminosity when matter is supplied to the system at a constant accretion 
rate. These are related to the periodic appearance and disappearance of the Eddington Capture Sphere. In the model 
the frequency is found to vary inversely with the luminosity. If the input accretion rate varies (strictly) periodically, 
the luminosity variation is quasi-periodic, and the quality factor is inversely proportional to the relative amplitude of 
mass accretion fluctuations, with its largest value Q ~ 1/(10 \8M/M\) attained in oscillations at about 1 to 2 kHz 
frequencies for a 2Mq star. 
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1. Introduction 



Abramowicz, Ellis & Lanza (1990) (hereafter AEL) argued 
that the luminosity of a relativistic star that is accret- 
ing at a super-Eddington rate, should periodically change. 
They have shown that in the combined gravitational and 
radiation fields of a spherical, compact star, radially mov- 
ing test particles are captured by a sphere on which the 
gravitational and radiative forces balance. This is because 
in Einstein's general relativity the radiative force dimin- 
ishes more strongly with the distance than the gravita- 
tional force, and radiation may be super-Eddington close 
to the star, but sub-Eddington further away, reaching 
the Eddington value at the Eddington Capture Sphere 
(ECS), whose radius is given by a simple expression in the 
Schwarzschild coordinates ( |Phinney 1987[ ), 
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Here R is the radius of the star, Rq = GM/c 2 its gravita- 
tional radius, and L, which is assumed to satisfy Eq. (|2|). 
is the stellar lum inos ity at its surface. [Bini et al. (2009)[ 



Oh et al (2011)[ and |Stahl et al. (2012)| have shown that 
the ECS captures particles from a wide class of non-radial 
orbits as welfl, so the following discussion is not restricted 



1 As shown in detail by Oh et al. (2011) azimuthal radiation 
drag efficiently removes particle's angular momentum, typically 
making motions in the combined gravitational and radiation 
fields asymptotically radial. 



to the case of radial accretion. However, it is important to 
keep in mind that the balance of forces necessary for the 
existence of the ECS requires a large, radial radiative flux. 
In non-spherical accretion this can only be realized in the 
optically thin regime. 

In this paper we present a simple model in which the 
idea of periodic luminosity changes is realized. AEL as- 
sumed that the radiation power was provided by the kinetic 
energy of the accreted particles, rele ased when they hit 



the surface of the star. As stressed by Stahl et al. (2012) 
particles settle on the ECS rather gently, so the ECS it 
self is not particularly luminous. If all particles are cap- 
tured at the Eddington sphere, they do not reach the 
surface of the star, and the stellar accretion luminosity 
goes to zero. When it does, and actually already when 
L/L Edd < (1 - 2i? G /i?)- 1 / 2 , the ECS disappears, accre- 
tion is resumed by the star and eventually the luminosity 
may build up to its former value, so the ECS will reappear 
and accretion will stop again. Thus, the accretion process 
may be quasipcriodic, alternating between states of high 
luminosity and no luminosity. 



2. The model 

The toy model considered here assumes a steadjO supply 
of optically thin flui d at some distance above the stellar 
surface. As shown by Stahl et al. (2012)] unless its velocity 
is extraordinarily large the fluid will settle on the ECS, 



The assumption of a constant M will be relaxed in Section[4] 
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regardless of its point of origin, if only 



(1 - 2R G /R)-^ 2 < L/L Edd < (1 - 2R G /R)-\ (2) 



In our model the stellar luminosity will be assumed to 
be either zero, or to have a certain definite value, L = Lq, 
satisfying Eq. ([2]). In this sense the model states are binary 
(on-off). This property implies that the ECS is located at a 
specific radius, Tecs — r o > R, whenever present (Eq. [Tj). 
The model is described by three equations: 



^ecs(*) 
L.(t) 

M„(t) 



fj,L*(t - 5t[), 

r]M*(t -St r ), 

M if ^Ecs(i - SU) < fj,L , 
if L E cs(t -SU) > LiL , 



(3) 
(4) 

(5) 



where Legs is the stellar luminosity at the ECS, L* is the 
luminosity at the stellar surface, M* is the accretion rate 
at the stellar surface, /i is a rcdshift factor between R and 
ro, r) is the conversion factor between accretion rate and 
stellar luminosity 

The simplest possibility, which we will ignore, is that 
r)Mo < ^Edd(l — 2Rg/R)~ 1 ^ 2 , so the critical luminosity for 
establishing the Eddington capture sphere is never reached, 
and the equations describe a steady state solution. We will 
assume instead that Lq = TjMo, i.e., the stellar luminosity 
in the "on" state has the necessary value to establish an 
Eddington Capture Sphere at the radius r . This assump- 
tion allows non-trivial time behaviour of the model system. 

There are three timescales that determine the behaviour 
of the model, the reaction time St r , which is the timescale 
for converting accreting matter into radiation, the light 
travel time Sti, which is the travel time of light from the 
stellar surface to the Eddington Capture Sphere, and the 
infall time Sti, which is the infall time of matter from the 
ECS to the stellar surface in the absence of radiation. Now, 
Sti and Sti can be determined from Rq,R and ro. So, for a 
given star these two timescales are a function of the stellar 
luminosity alone, i.e., of r/Mo- The reaction time introduces 
a phase shift between accretion and radiation at the stellar 
surface. In numerical examples we will adopt two extreme 
values St r = 0, or 5t r = 0.1ms, the latter being the esti- 
mated cooling time of a clump of fluid that falls on a neu- 
tron star surface (Kluznia k, Michelson, fc Wagoner, 1990[ ). 

In reality, the time delays would not be as sharp as we 
have assumed them to be. For instance, even at time Sti 
after the luminosity is turned off, particles outside the ECS 
suffer from radiation drag for an additional interval of time 
(r — Tq)/c, where r is the position of the particle at time 
Sti + (' — r o)/c after the radiation is turned off. Conversely, 
after the luminosity is turned on again, all particles between 
the ECS and the stellar surface will feel the full impact of 
radiation pressure after a delay of only (r — R)/c < Sti, 
where again, r is the position of the particle when the radi- 
ation front passes it, (r — R)/c after the radiation is turned 
on. In the toy model, we neglect such effects entirely. 

Combining Eqs. ©-I® wc find that M*(i) is determined 
by M«(t - T), with T = St, + St r + Stf. 
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Fig. 1. Stellar accretion rate and luminosities at the stellar 
surface and at the Eddington Capture Sphere, according to 
the model. 
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Fig. 2. Stellar accretion rate and the luminosity at the 
ECS, for another choice of initial conditions. 



M«.(f) = M Q(jiLo-LEGs(t-5U)) 
= M G{L - L*(t - St t - Sti)) 

= M e (m - M*{t - su - sti - St r ) 

= Mo-M^t-T), 

where is the Heaviside step function. Thus, 

{M*(i -T) = Mo} => (m*(*) = o} , 

{M*(f — T) = j => {M*(t) = M } , 

and clearly, the model system shows periodic behaviour 
with the period 2T: 



M*(t) = AU{t-2T) 



(6) 



3. Results of the model 



A typical behaviour of the model is shown in Fig. Q] which 
shows a solution of Eqs. (OJ)- (J5j) with Sti ■ Sti ■ St r = 8 : 
17 : 10. The black solid line traces M*. In the figure, we 
identify the time b — a = e — d = St r as the reaction time, 



2 



Maciej Wielgus et al.: Oscillations of the Eddington Capture Sphere 




Fig. 3. The semi-period T of the oscillation in geometri- 
cal units as a function of peak luminosity at infinity in 
Eddington units, when kinetic energy is instantaneously 
converted to luminosity (5t r = 0). The corresponding fre- 
quency / = 1/ (2T) can be read off in kHz from the right 
vertical axis. The curves are labeled with the stellar ra- 
dius (3, ..,10) in units of Rq- Filled squares indicate a value 
of the slope dlogT/dlogLoo = 50, filled circles a value of 
100, and the crosses the minimum value of the logarithmic 
derivative for each curve. 



c — b = h — e = Sti as the light travel time, and d — c = 
j — h = Sti as the infall time. The luminosities at the stellar 
surface and at the ECS are shifted relative to the accretion 
rate by St r and St r + Sti , respectively. We note also that 
d — a = j — d = h — c = k — e = e — b = St r + Sti + Sti = T, 
and the pulses are as long as the periods without activity. 
Thus the period of the oscillation is P = 2T, as expected 
from Eq ((6]). 

In constructing Fig. [T] it was assumed that matter first 
arrives at the stellar surface at t — and that there is a 
continual inflow of matter to the system at the rate Mq 
for all t > 0. Different initial conditions can lead to a 
more complicated pulse shape of M*(i) for t in the in- 
tervals (2nT,2nT + T), with its complement M„(t) = 
M - M* (t - T) in the intervals (2nT + T, 2nT + 2T). Here, 
and elsewhere, n = 0, 1, 2, 3.... The state of accretion (Mo or 
0) at time t has no influence on the future value of M» until 
the time t + T, so one can assume that at any instant in the 
initial interval t £ [0,T) the accretion rate has any of the 
two values, or Mq, i.e., in this interval M*{t) = Mog(t), 
where g(t) is an arbitrary binary function, mapping the 
initial time interval into on-off states, g : [0,T) — > {0,1}. 
Fig.[Uprovidcs an example of such behaviour. Thus, in prin- 
ciple, the harmonic content of the signal may be quite rich, 
although the fundamental is still at / = 1/P = 1/(2T). 

As can be seen from Eqs. (H])-© the model takes no 
account of accumulation of matter on the ECS, that is to 
say, here we ignore the arrival at t = 2(n + 1)T (and at 
other moments as well in the case of Fig. [5J , of the addi- 
tional matter that had accumulated on the ECS. In reality, 
this matter would produce a brief flash of radiation of en- 
ergy ?7ecsM3AT, where t?ecs is the conversion efficiency 
to radiation of the kinetic energy of matter falling from 
the ECS, and AT is the accumulation time (AT = T in 
Fig. [J). If in the steady accretion phase matter is falling in 
from r >> rg, one would expect ?7ecs << V- 



Fig. 4. The frequency of the oscillations as a function of 
peak luminosity at infinity in Eddington units, when there 
is a delay in converting kinetic energy to luminosity St r = 
0.1ms (see Kluzniak, Michelson, fc Wagoner, 1990[ ). The 



stellar mass is assumed to be M = 2M Q . The curves are la- 
beled with the stellar radius (3,. .,10) in units of Rq- Filled 
squares indicate a value of the slope d log T/d log = 30, 
filled circles a value of 50, and the crosses the minimum 
value of the slope for each curve (boxed values). 



In a physically realistic situation the infall time is the 
dominant time scale, Sti ^> Sti. The radius of the Eddington 
Capture Sphere, and hence also the infall time, increases 
rapidly with L^, and without bounds as — > L^dd- 
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Here, L M is the redshifted luminosity at infinity. Hence, 
except for the lowest values of the luminosity parameter 
Lq, the infall time Sti is the dominant time scale, and the 
frequency of oscillations / = 1/(2T) varies inversely with 
the luminosity. Neglecting St r , we can compute the semi- 
period T, and the frequency of oscillation as a function of 
the stellar radius R and of the luminosity L M alone — T 
and / will, respectively, scale directly or inversely with M. 
These are shown in Fig. [3] for various stellar radii. Fig. 2] 
shows the frequency as a function of luminosity, for a star 
with M = 2M Q , when St r = 0.1 ms. 

4. Conclusions and discussion 

We have shown that supplying mass to the vicinity of a 
compact star, continuously and at a constant accretion 
rate, can lead to periodic top-hat variations of luminosity, if 
only the mass accretion rate corresponds to a mildly super- 
Eddington luminosity at the stellar surface. This oscillatory 
behaviour of luminosity is related to the phenomenon of the 
Eddington Capture Sphere, which is a consequence of the 
interplay of radiation drag and general relativity. 

If such an oscillation occurs in the real world, for ex- 
ample in the Z sources, where rapid variations of the in- 
ferred inner radius of the accretion disk have been reported 
(ILin, Remillard, & Homan, 20091) , the actual mechanism is 



likely to be more complex than the strictly periodic oscil- 
lation in the toy model considered here. E.g., the accretion 
rate is not likely to be constant (contrary to our assumption 
of simple on-off behaviour at the stellar surface) . 



3 



Maciej Wielgus et al.: Oscillations of the Eddington Capture Sphere 



As a minor extension of the model, consider an accretion 
rate that is varying on a timescale comparable to the period 
of the ECS oscillator. Even strictly periodic variation of Mo 
will lead to a decoherence of the ECS oscillation, because 
the position tq of the ECS and (hence) the infall time from 
the ECS, and (hence) the oscillator frequency, are strongly 
varying functions of the luminosity, Lq = rjMo- Indeed, 
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Thus, the Q factor of the oscillation is 



Q 
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(8) 



(9) 



We have indicated some values of the logarithmic derivative 
d In T/d In L m in Figs. [31 |U Comparing the figures, we see 
that delayed emission (reaction time St T > 0) at the stellar 
surface acts as a low-pass filter, cutting out the highest 
frequencies. At the same time it improves the quality factor 
of the oscillator. 

Accreting neutrons stars in low mass X-ray binaries are 
expected to have a radius in the range of 4 to IGM/c 2 
( |Arnett fc Bowers, 19771 |Kluzniak fc Wa goner, 19856, and 
a mass of about two solar masses. It is seen from Fig.|4]that 
the minimum values of the logarithmic derivative attained 
for radii of 4, 5, 6, IGMjc 2 are between 5 and 15. Hence, 
the maximum expected value of the oscillator quality factor 
would be Q sa l/(10|£Mn/Mn|), occurring for frequencies 
between about 1 and 2 kHz, as can be seen from Fig. [4] 
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